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Abstract 

General properties of the light and heavy quark propagators have been investigated 
in the context of string tension interaction. Confinement, chiral symmetry breaking, 
spectral properties of the propagator are analytically studied and numerically validated. 
We show that the propagator is analytic in the infrared region even for massless quarks 
with a non zero radius of convergence. Emergence of more than one mass scale is 
exemplified. Massless limit of the quark propagator does exhibit critical behaviour. 
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1 Introduction 



Confinement of quarks and their transient behaviour as dictated by the quark propagator is 
being investigated in a theory jTj which shows many features of QCD. We have considered 
colored quarks in QCD with renormalisable string tension in the large N expansion. We 
compute the behaviour of the quark propagator from first principles in this quantum field 
theory as a function of momentum, renormalised mass and coupling constant. 

The presence of string tension results in the absence of the pole in the quark propagator 
thus manifesting the non-existence of quarks as asymptotic states in the quantum field 
theory. In addition for any renormalised mass parameter m > 0, we find that the analytic 
structure of the propagator is non-trivial. Firstly for space like momenta (p 2 < 0) the quark 
propagator is analytic. For any renormalised mass there exists a Taylor series solution in p 2 
to the Schwinger-Dyson(SD) equation. This series is numerically estimated to converge for 
p 2 < fh 2 ^ where fh is called the threshold mass. In all, associated with a quark, there 
are three mass parameters namely, the renormalised mass m, the threshold mass fh and 
the asymptotic mass M which capture the behaviour of the quark propagator for space like 
region. These three masses are very distinct and are estimated from phenomenology. We also 
discuss how the so called constituent mass in various bound states is related to fh. We find 
for heavy quarks where the renormalised mass scale m 2 is larger than string tension a, there 
is a qualitative difference in the analytic properties which can be physically interpreted as, 
in space-time light quarks have zero transient time existence while heavy quarks are almost 
about to be free and hence have relatively long time transient existence. 

Massless quark theory does exhibit spontaneous chiral symmetry breaking in the pres- 
ence of string tension and the theory is self-consistently defined with out the need for ul- 
traviolet (UV) renormalisation unlike the massive case. Consequently the quark propagator 
function is not a continuos function of m at m = 0. The spin independent quark propagator, 
due to Goldstone phenomenon, is also the massless pion wave function with the asymptotic 
mass of the quark in spacelike region being the inverse of the effective size of the pion. The 
SD equations also admit a chiral symmetric solution with confinement, this is understood to 
be an unstable vaccum. We also show that photon like behaviour of the gluon cannot cause 
chiral symmetry breaking for any value of the strong coupling constant g in the rainbow 
graph approximation. 

Much of this analysis is done analytically by exploiting certain computational simplifi- 
cations that we state in the paper. This simplification procedure can be applied for any 
relativistic field theory. Numerical investigation was essential to establish the existence of 
solutions in the space like region and thus estimate the asymptotic masses of the quarks as 
a function of the parameters of the theory. A recent review on SD equations is given in [2|. 

2 Lagrangian with String Tension 

The Minkowski Lagrangian of the theory is given by [Tj 
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Z = J DADQDCDxD X exp(iS + i J j.(Q + A)) } 

So = J (^F 2 + C,(-D 2 )Q^ + X,(-D 2 ) X " ~ ~CJ) + / (qi^D»(A + Q)q - mqq) (1) 

where F = F£ v = d^A^ — d v A a ^ + f ahc A h A c v is the antisymmetric gauge field tensor, A ^ are 
the gauge fields, g is the gauge coupling constant, f abc are the structure constants of the 
corresponding non-abelian gauge group. and are bosonic vector fields and x an d 
X are Grassmann Lorentz vector fields. All of them are in adjoint representation of the 
gauge group and transform covariantly under local gauge transformations. The significance 
of these fields is discussed in [T]. D ab = D^A) = d^8 ab + f acb A c ^ is the covariant derivative 
operator, q and q are quark fields (fermions) in the fundamental representation of the gauge 
group with a mass m. As can be seen from the Lagrangian the quark fields have a qAq and 
a qQq interaction vertices, a is the string tension which is asymptotically free jT]. Gauge 
fixing for the gluons has to be done just as in QCD and in the Feynman gauge the A^ and 
the propagators are given by 

F {q z + ie) z * q 2 + xe 

We consider a systematic expansion in the various coupling constants aN,g 2 N and 1/N 
, where a is the string tension, g is the gauge coupling and iV is the number of colours, of 
which the leading infrared(IR) contributions come from the aN term alone. Consequently we 
perform a non-perturbative expansion in aN but perturbative expansion in g 2 N and 1/N [3J. 
This yields for the quark propagator a summation of the rainbow graphs with the < Q^Q^ > 
propagator alone . The SD equation can be considered with only this term alone and solved 
self-consistently. It turns out this equation has no UV infinite renormalisation. But we know 
that including the g 2 N term requires UV renormalisation. This raises the question whether 
non-perturbative IR physics is truly independent of the UV renormalisation procedure. To 
address this question in our analysis we include in the SD equation all the leading contri- 
butions in the various asymptotic IR and UV regions. The IR region is dominated by the 
< QuQ b > propagator whereas the UV region is dominated by the < A° il A b v > propagator 
which gives the third term in the SD equation. It should be mentioned that to order g 2 N in 
addition to the above term there are some more graphs but these are less important in the 
UV region and are not expected to make a significant contribution, for simplicity we neglect 
them. In this work we do not use running a, g or m in our calculations. The SD equation 
for the quark propagator with both the < Q^Ql > and < A^A b u > propagator (Fig[T]) is, 

1 r rt A k aN n 2 N 

S - Hp) = _ (7 . p - m) + / _ ( __ 7 „ S(p - k)Y - jjfcyM - k hn (3) 
where S(p) is the full quark propagator with momentum and mass m. 
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Figure 1: Schwinger Dyson Equation for Quark propagator 



The dashed lines in the figure corresponds to < Q a u Q b u > propagator and the springy 



line to the gluon.The solid line with the blob is the full quark propagator and the solid line 
without the blob is the bare quark propagator. Spin decomposing the quark propagator S(p) 
as S(p) = i(A{j> 2 )^ .p + B{p 2 )) and substituting it in Q, we can express the SD equation as 
a coupled nonlinear equation in terms of two scalar functions A(p 2 ) and B(p 2 ). 

3 Causal Representation 

The solution space of SD nonlinear equations in general can be very large to address. Rela- 
tivistic Quantum Field Theory does impose certain qualifications and hence only those alone 
are of interest to us. Generally these restrictions are encoded in the Kallen-Lehmann repre- 
sentation. The physics consideration that go into this are that positive energies propagate 
forward in time and negative energies backward in time (as imposed by the Feynman ie pre- 
scription) and all intermediate states have a positive norm. While these are essential for any 
physical propagators, it remains questionable whether unobservables should also have these 
properties. In particular the p propagators necessarily generate negative norm intermediate 
states in perturbation theory 0]. Therefore we disregard the second assumption and make 
the minimal assumptions necessary namely Lorentz covariance and Feynman ie prescription 
yielding the causal representation. 



A(p 2 ) = r da 2 A{a) , (4) 
J-oo p 2 — a + ie 

/ 9n f°° , B(a) , , 

B(p 2 ) = / da 1 L . 5 
J-oo p 2 — a + it 

where A(a) and B(a) are functions that are not required to be positive. Infact the above are 
the Cauchy-Riemann integral representations for analytic functions. All we are demanding 
is that the quark propagator that we are interested in should be an analytic function of 
momentum. From the analyticity properties one can infer that the imaginary parts of A(p 2 ) 
and B{p 2 ) are the causal weight functions A(a) and B(a) themselves. In eq.(dJ) and eq.© 
out of generality we consider the entire range of a. However as it turns out in the succeeding 
analysis A(a) and B(a) are non-vanishing only for a > 0. 
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4 Regularisations 

The < Q^Ql > propagator behaviour l/(k 2 + ie) 2 is highly singular in the IR regime. In 
3 + 1 dimensions the +ie prescription is not sufficient to regulate it. Even after Wick rotation 
it has a logarithmic singularity which makes this propagator ill-defined. Consequently the 
theory demands a regularisation of this IR divergence. There are many equivalent ways 
of regularising this but we will mention the one which is convenient to implement in our 
relativistic quantum field theory. In all loop calculations we are convoluting this singular 
propagator with another quantity which has an integral representation which in its generic 
form is ^ k _ p ^2_ a+ie an d higher powers thereof. We illustrate the IR regularisation prescription 
by explicitly showing the steps in carrying out the B(p 2 ) integration. Consider the following 
integral that occurs in the SD equation, 

d 4 k B((p-k) 2 ) 



i(2vr) 4 k A 



(6) 



In eq.flBJ, we substitute the spectral representation eq.(jEJ) for B((p — k) 2 and using Feynman 
parametrisation, we do the 1-loop momentum integral explicitly by making the standard 
Wick rotation. 

30 r d 4 k 1 B(a) r°° da r 1 dx xB(a) 

oo J z(2tt) 4 {k 2 + ie) 2 {k - p) 2 - a + le = J-oo J^f J (1 - x) [xp 2 - a + ief ' 

<lv xB(xp 2 ) (8) 



(47r) 2 Jo (1- x 

In the above expression we have done the a integration formally to get back B(p 2 ). As is 
evident in eq.(jHJ) the IR singularity manifests itself as the end-point singularity at x — 1 in 
the Feynman parameter x. The divergence is regulated by explicitly subtracting the x = 1 
singularity in the integrand. 

dx dx 

/ n \x b (xp 2 ) - / n -AxB{xp 2 ) - B{p 2 )) (9) 

Jo (1 — x) Jo (1 — X) 



f 1 d 
/ dx ln(l — x) — (xB(xp 2 )) 
Jo dx 



This defines the regularisation which can be unambiguously implemented in the theory. The 
same procedure follows through for the integration involving A(p 2 ). 

The gluon term in the SD equation is UV divergent which needs to be regularised. The 
standard dimensional regularisation and subsequent multiplicative renormalisation of the 
quantum fields and parameters can be implemented. We find the same can also be obtained 
in a uniform way by implementing the following representation for the gluon propagator 

I poo l 

l-r = - dP (10) 
p l + ie Jo (p 2 — p + it) 
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Using the causal representation we can do the loop momentum integration. The UV divergent 
terms are f °° d/3A(j3) and J °° df3B(f3) which can be removed by renormalising the wave 
function and the quark mass parameter as in ordinary perturbation theory. The finite SD 
equation for renormalised mass m are, 



A(p 2 ) r 1 x 2 A(xp 2 ) - A(p 2 ) 



? A( ?y - B{p 2 ) 2 l + 2d L dx * Z-x) - ap2 1 dxA{xp2){x ' 2 - 1] 
pM(p |^W = m + 4g /.' ^fl^ - [ W>(* - 1) (») 



where we have defined a = and a = . In the rest of the paper for convenience we 
will work in mass units of a = 1. 



5 Confinement and Asymptotics 

All states which exist for asymptotic times are physical and they manifest as poles in the 
corresponding Green's functions. Physical quarks should be realised as poles in the quark 
propagators. Looking for poles in eq. (jllj) . if A(p 2 ) has a pole for some value of p 2 , the 
left hand side of eq.([ll|) vanishes and the right hand side using the last representation in 
eq.(fTU)l. shows an edge singularity divergence. This contradiction in eq.()lip implies absence 
of pole like solutions. Similar is the case for B(p 2 ). Finally if both A(p 2 ) and B(p 2 ) have 
poles at the same value of p 2 then by taking the ratio of eq. pijl and eq. ()12)) we note that 
again the equation cannot be matched. These observations make us conclude that A(p 2 )&nd 
B(p 2 ) cannot have poles, furthermore we conclude by refining this observation that A(p 2 ) 
and B(p 2 ) cannot even have any divergent behaviour for any value of p 2 . 

To analyse the large space like p 2 (p 2 — > — oo) behaviour of A(p 2 ) and B(p 2 ) in eq.tjllj) 
and eq.(jI2J), we rescale the integrand x — > and obtain the approximate equation, 

1 rP 2 ^A(x) - A(p 2 ) rp 2 x 2 

1 + 2-/ dx p )' - — -a / dxA(x)(--l) (13) 
p 2 Jo (1-^) Jy p A 

1 rP 2 ^B(x)-B(p 2 ) rp 2 x 

m + 4-/ dx p / yFJ -4a dxB(x)(- - 1) (14) 
P 2 Jo (1 -3) Jo ,y p 2 
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p 2 A(p 2 ) 
B{p 2 ) 

2\2 



p 2 A(p 



We further split the integral range to (0, 1) and (l,p 2 ). It is easy to see that the latter dom- 
inates and we can determine the asymptotic behaviour of A(p 2 ) and B(p 2 ) self-consistently 
in the UV regime. It is interesting to note that the UV behaviour of A(p 2 ) and B(p 2 ) is not 
influenced by the IR behaviour at all as expected in perturbation theory. 
For p 2 — > —00 

A{p 2 ) ~ 1(1, 1 ) (15) 

P 2 J2a\n(-p 2 ) 
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B(p 2 ) ~ \(m, ln(-p 2 )) (16) 



In the above equations the first term in the bracket refers to the asymptotic behaviour in 
the absence of the a term while the second term refers to the asymptotic behaviour in the 
presence of the a term. It has to be noted that except in the chiral limit, upto logarithms the 
large momentum behaviour of A(p 2 ) and B(p 2 ) with and without the gluons are the same. 

Next we look at the small momentum (p 2 — > 0) limit. From the integral equation the 
small p 2 behaviour of A(p 2 ) and B(p 2 ) are also self-consistently determined without any 
contribution from the UV region. Infact at p 2 = 0, A(0) and B(0) are determined to be 



B(0) = ± V f + 16 (17) 

.4(0) = 8fl '°» ■ 
v ; m-B(0) 

Near p 2 = we can have series expansion solution A(p 2 ) = Y^=o a nP 2n and B(p 2 ) = 
En=o b nP 2n for eq.(HU) and eq.jEJ. Substituting the series and equating the coefficients of 
like powers of p 2n , we get the following recursion relations for the coeffecients b n and a n , 



a, 



bn ~\~ 2 X/m=0 (.^^"mbn—mfn—m < ^(20/ m 6 n _ m _i/l n _ m _i d n — m —ib m g n — m —l)) 

m - Ab f + 2b f n +i 

2 X/m=l Q"n—mbmfn—m+l 

m - 4b f + 2b f n +i 



7 0>n—\ 2 5Z m= o(o , m'^n— m— l/n— m "I - 2oft m ft n _ m _i/l n — m _i) -|- 2tt X/m=0 ^"m^n—m—29n—in—2 

m - 4b fo - 4ft / n 

4V n_1 ft ft f 

_|_ ^ Z-im=l u m u n—mJ n—m (19) 

m - 4&0/0 - 4& /„ 

where, 

i 1 1 1 

/i = §t+T ;/ll= (* + l)(* + 2) ; ^ = (i + l)(j + 3) 

The radius of convergence of the series solution is difficult to ascertain analytically. In- 
vestigating it numerically, we find that for small masses m, the series converges in a small 
interval around p 2 = 0. As we increase the mass m, the interval of convergence also increases. 
We attribute the finite radius of convergence to the onset of the imaginary part in the SD 
equation at some value of the p 2 > and this value is proportional to the mass m. The 
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onset of this behaviour is perhaps due to a branch cut and consequently no Taylor expansion 
can converge. Numerical solution of the integral equations eq. (jll|) and eq. (|12|) support this 
assumption. The onset of non-analyticity in A(p 2 ) and B(p 2 ) in general can be different but 
we find that for small quark mass m they are about the same. We cannot ascertain precisely 
the nature of non-analyticity but it can be A{p 2 ) ~ const + w^pz^ or even softer. This 
follows because if we assume that A(p 2 ) or B(p 2 ) is singular at some value of p 2 , the r.h.s of 
eq. ljll)) and eq. (jl2)) will yield even more singular contribution for the same value of p 2 thus 
mismatching the equation. Similar is the case for B(p 2 ). The Taylor expansion suggests that 
for small quark masses m, for both A(p 2 ) and B(p 2 ) functions the non-analyticity threshold 
is at m 2 and is numerically estimated to be 

fh 2 ~ .02 + m * const for m < 1 (20) 
fh 2 ~ m 2 for m > 1 

For large quark masses we estimate m 2 by doing the naive large m or equivalently the small 
a and a perturbation theory. In this perturbation theory where a is also small , it is easy to 
recognise that A(p 2 ) and B(p 2 ) may not have strict poles but poles softened by logarithms 
such as ■ Hence there is a net singularity softer than a pole. This is a remarkably 

different behaviour as compared to that of light quarks. 

In eq. (jl7p . if we pick B(0) negative solution, we are able to establish analytically that the 
imaginary part may exist only for p 2 > and both A(p 2 ) and B(p 2 ) are negative for space 
like p 2 . From these two observations it is self-consistent to infer that that the imaginary 
part of A{p 2 ) and B{p 2 ) are also positive. Consequently we can conclude that the quark 
propagator in our theory is consistent with the standard Kallen-Lehmann representation. 



6 Numerical Solutions 

The numerical solutions for A(p 2 ) and B(p 2 ) in the space-like (p 2 < 0) region can be obtained 
by iteration. Solutions exist with or without the gluon term. We ignore the gluon term for 
simplicity . The functions A(p 2 ) and B(p 2 ) are recast in terms of two new functions A(p 2 ) 
and B(p 2 ) through the following transformation 



■U/r 



A(j) ^ p 2 A(p 2 ) 2 - B(p 2 ) 2 

U(n 2 \ = (91 1 

KP ' p 2 A{p 2 ) 2 - B{p 2 ) 2 y ' 

This is substituted in eq. ljlip and solved for A(p 2 ) and B(p 2 ) by giving a seed which is 
A(0) and -B(O). Such an equation is numerically convenient to iterate. A(p 2 ) and B(p 2 ) 
can be obtained by solving eq. (}2"Tj) in terms of A(p 2 ) and B(p 2 ). The nature of the integral 
equation !) 11)1 is such that to solve it at a particular value of momentum, one has to know the 
functions at all values of momentum below that particular value all the way upto p 2 = 0. 
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The momentum p 2 is discretised and the integration over x is done numerically. We use an 
interpolation formula to fit the iterate and do the x integration. For values of m > 0.5 the 
iteration method works very well. But for small values of the quark mass m, due to the 
onset of the imaginary part very close to p 2 = 0, one has to be very careful. The solutions 
for very light masses have been determined with about 5% inaccuracy. 

The solution for light quark masses are shown in Fig[2K] and Fig(2b]- The string tension 
a sets the scale for small p 2 and the dominant contribution in this region comes from this 
term. For p 2 > 1 the momentum p 2 sets the scale. The large momentum behaviour given 
by eq. (jl5|) and eq. (jl6p sets in at around p 2 ~ 3. For small values of m, A(p 2 ) has a minima. 
For masses m > 0.82 this minimum does not exist in the spacelike region. 
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Figure 2: Numerical solutions for (a) B(p 2 ) and (b) A(p 2 ) for light masses 



The solutions for heavy quark masses, are shown in Fig[3k] and Fig|2t>]. Unlike in the 
case of light quarks, the mass m dominates over the string tension and sets the scale for small 
value of momentum. For large momentum like in the light quark case, its the momentum 
that sets the scale. 

It can be noted that all the solutions for A(p 2 ) and B(p 2 ) in the large momentum limit 
goes over smoothly to the asymptotic expressions in (JT5j) and (jl6|) . Upto to logarithmic 
corrections, A(p 2 ) and B(p 2 ) have essentially the canonical p 2^ M 2 behaviour in the space- 
like region for all masses m where M is called the asymptotic mass. A useful parametrisation 
for A(p 2 ) and B(p 2 ) is, 



A{p 2 



-A(0)M 2 



M 2 



B(p 2 ) ~ 
{P ' p 2 -M 2 



(22) 
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Figure 3: Numerical solutions for (a) B(p 2 ) and (b) A(p 2 ) for Heavy masses 



The "asymptotic masses" M 2 and M 2 depend on the mass m. A convenient way of 
determining M 2 and is by defining them to be that value of p 2 where A(p 2 ) and B(p 2 ) 
are reduced by half from their values at p 2 = 0. Doing so we can infer the asymptotic masses 
M a and Mb from the numerical data. In Fig^j we have plotted M a and Mj, versus m and 
for comparison we have also plotted M = -^jjjy which is 



Mo = 7m + + 16 (23) 
8 

The asymptotic mass does increase almost linearly with m, all through from light to heavy 
quarks. There are two asymptotic masses for any quark, namely M a , the spin-dependent 
asymptotic mass and Mb the spin-independent asymptotic mass and the former is signifi- 
cantly heavier than the latter for light quarks. 

We can extend the numerical ansatz to small p 2 > as well. Beyond which as we 
understand from the Taylor expansion the numerical solutions cannot converge . At this 
point we have to work with the real and imaginary parts of A(p 2 ) and B(p 2 ) which is 
numerically cumbersome to handle. The asymptotic mass behaviour of the quark propagator 
given in eq. (|22j) is to be taken for spacelike p 2 only. For timelike p 2 > this is qualitatively 
unacceptable due to confinement. 

The above procedure of determining numerical solutions can be done for B(0) positive 
also. We find that the numerical solution is not so stable. This is partly because in the the 
r.h.s of eq.(|12|). we can see that the integral contribution is of opposite sign to that of the 
constant term m. Therefore at some point B(p 2 ) will vanish and change sign. The existence 
of this solution can be doubted numerically but it cannot be ascertained. In any case this 
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Figure 4: Asymptotic mass M a , Mf,,M for A(p 2 ) and B(p 2 ) v/s renormalised mass 



solution is not of physical interest as it disobeys Kallen-Lehmann positivity condition. From 
eq.(|l*|) and eq.l j^j l. knowing that A(a) and B(a) are non vanishing for p 2 > m 2 > 0, A(p 2 ) 
and B(p 2 ) will never vanish in the space like region if A(a) and B(a) are positive. 



7 Chiral Limit 



We look at the theory when the bare mass of the quark vanishes. Then the theory has a 
global continuos symmetry namely chiral symmetry. In our large N and small g 2 N expansion, 
doing just the wave function renormalisation, we obtain the SD equations where we have 
eq.(JTTJ) and 



The last term has a constant perhaps finite or infinite. If it is infinite, then renormalisation 
is necessarily demanded. If it is finite renormalisation is an option. In either cases if we do 
any renormalisation, finite or infinite, it is evident that this implies we have formally included 
a bare mass term and then subsequently arranged the renormalised mass to vanish. This 
procedure cannot respect chiral symmetry. This is easy to see by noting that the massless 
Goldstone meson wave function satisfies the following Bethe-Salpeter equation jH], 

0(p») = (p 2 A(p 2 ) 2 - B(p 2 ) 2 )(4a f 1 dx X<P{x fP - f (p2) + 4a [' dx d^(P)) (25) 

JO {1 — X) Jo J~oo 

In a renormalisable theory, <f)(p 2 ) does not undergo any other renormalisation. Conse- 
quently in eq.(j2HJ the last term cannot be altered. It is clear that the solution to eq.(j25J) 
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is </>(p 2 ) = B(p 2 ) provided we do not make any finite or infinite renormalisation of eq.([24|). 
4>{p 2 ) = B{p 2 ) is precisely the consequence of the Goldstone theorem. This property of SD 
and Bethe-Salpeter equations, first realised by Mandelstam is true in theories with spon- 
taneous symmetry breaking in the rainbow graph approximation. To reiterate, any mass 
renormalisation scheme if adopted such that the constant term in the last term in eq. ()24|) is 
effectively absent then chiral symmetry is explicitly broken [HI U\ ■ 

Alternatively we can investigate if eq. (J24)) as such has any solutions. The crucial property 
for the existence of the solution is its UV behaviour. Using exactly the same approximation 
as in eq. ()13|) and eq. (fi4"j) we find that behaviour of A(p 2 ) for large spacelike momentum is 
unchanged from eq.(fl3|). while B{p 2 ) behaves as 



B(p 2 ) ~ C ° nst ( 1 1 \ ( 2 Q) 

{p 2 ) 2 p 2 ' (ln(— p 2 )) 3 

The UV behaviour of B(p 2 ) is indeed more convergent than the canonical behaviour. Con- 
sequently the last term in eq. (|2"4"|) is finite. Thus our system is self consistently determined 
with only wave function renormalisation and no mass renormalisation. 

In the small p 2 regime not much of a difference occurs and indeed -B(O) = ±| ( neglecting 
the a term) as expected from eq. (JTTJ). Asymptotic behaviour of B(p 2 ) as in eq.(|2^j) is 
determined upto a constant and this can have either sign. Now it is possible to interpolate 
both the solutions to their corresponding asymptotic behaviour without B{p 2 ) vanishing. 
These solutions are chiral conjugate. This is because in eq. (|24|) there is an ambiguity in 
the sign of B(p 2 ). Both B(p 2 ) and —B(p 2 ) are possible. Thus picking one solution breaks 
the chiral symmetry spontaneously. Note that A(p 2 ) remains the same irrespective of the 
choice of B(p 2 ), because by definition A(p 2 ) is chirally symmetric. Both A(p 2 ) and B(p 2 ) 
have analytic expansions around p 2 = as determined from the series solution, eq. ()18|) and 
eq. ()19)) with a = 0. The radius of convergence of this series is numerically inferred to be 
nonvanishing as shown in eq. ()19)) for a = 0. 

Since the asymptotic behaviour of B(p 2 ) is very different from the canonical, we should 
expect the asymptotic mass definition to be modified. Namely for space like p 2 



B(0)(M 2 ) 2 -M 2 1 , (c?7] 

KP ' ~ (p 2 - M 2 ) 2K (p 2 - M 2 Y (ln(-p 2 )) 3; K ' 

Numerically the asymptotic mass can be estimated and it is found that M 2 ~ (0.7 — 0.9). 
The qualitatively different behaviour of B(p 2 ) in the exact chiral limit cannot be obtained 
continuosly by taking the renormalised mass m — > 0. Hence this theory exhibits critical 
behaviour at m = 0. This as we see is a consequence coming totally due to the absence of 
mass renormalisation in contrast to the massive case. The chiral limit exhibits spontaneous 
symmetry breaking but also shows that the spin independent propagation of a quark is non- 
analytic as a function of the quark mass m at m = 0. Inspite of this in the IR region as seen 
from the Taylor series we do not see any signature of the critical behaviour. 

SD equations (jllj) and (J24"| also have chiral symmetric solutions wherein B(p 2 ) vanishes 
for all p 2 and A(p 2 ) is non-vanishing. We can show there exists a series solution for A(p 2 ) = 
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J2^=o Q"a{~ P 2 )^ for small p 2 and the asymptotic behaviour remains unchanged as in 



(|15l) . This solution preserves chiral symmetry. Note that this is also a confining solution as 
there is no pole for the quark propagator. Indeed the propagator is finite at p 2 = and has 
only a square root branch cut starting from the origin i.e., the threshold mass vanishes. SD 
equations can have many solutions. Each of them actually correspond to different possible 
stable or unstable vacua. The criterion to pick the minimum vaccum configuration cannot 
be inferred from the SD equations alone. This can be carried out by summing the vaccum 
graphs which in turn depend on the propagator functions. We have not attempted this 
analysis here. By our expectations either we have one solution to the SD equation in which 
case that will correspond to the vaccum or there are three (odd) solutions of which two are 
related by chiral symmetry and hence to be degenerate stable vacua and the third has to be 
unstable. Any choice of the stable vacua therefore will yield symmetry breakdown and the 
chirally symmetric choice is unstable. 

We now address symmetry breakdown between two stable vacua in the absence of string 
tension. In eq. ()24|) we put o = 0. This is a truncation of QCD in which we are keeping 
the rainbow graphs due to gluons alone. This is not a controlled approximation but it is 
interesting to know whether less singular interactions with sufficiently large coupling constant 
a can cause spontaneous symmetry breakdown. The answer as we prove now is that it cannot 
cause symmetry breakdown. The asymptotic UV behaviour remains the same as before. In 
the absence of the a term we do expect the quark propagator to obey the Kallen-Lehmann 
representation, consequently for space like p 2 , A(p 2 ) and B(p 2 ) cannot have any other zeroes 
other than at asymptotic limit. Furthermore we also know if -B(O) is positive (negative) 
and it remains positive (negative) in the entire space like regime, due to Kallen-Lehmann 
representation. The consequence in eq. ()24|) at p 2 = is that the r.h.s is negative (positive) but 
the l.h.s is positive(negative) showing an inconsistency. Thus the only consistent alternative 
is -B(O) =0 and hence no chiral symmetry breakdown. 

8 Discussion 

We speculate some consequences of our semi-analytic analysis supported by numerical es- 
timates for strong interaction phenomenology. In [8J by fitting the vector meson mass we 
estimate m ~ .5Gev. This threshold mass is not a physical mass accessible at asymptotic 
times, however in most dynamics of quarks it still should give similar effects in the transient 
existence. For example in stable bound states quarks would have to have energy less than 
m. Consequently an electromagnetic probe will find the quark degrees of freedom to have a 
"constituent" mass depending upon the bound state but necessarily less than the threshold 
mass. Similarly a high energy quark will radiate and lose energy but it will cease to radiate 
when it reaches the energy of the threshold mass. 

From jH] the light quark (u and d) renormalised mass is estimated from the pion mass 
to be about 6Mev . This is consistent with other phenomeno logical considerations [5] and is 
substantially different from m. In this theory string tension runs and it decreases at high 
energies. We define heavy quark mass in comparison to relevant string tension. If m 2 > a, 
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we find that all notions of masses m, rh and M are the same. Furthermore heavy quarks 
in some sense are not severely effected by string tension. Apart from confinement, it is 
suggested that heavy quarks have a longer transient existence in contrast to light quarks. 

String tension, naively can cause non-unitary behaviour in the theory. We find that with 
in this limited exercise of quark propagator and their bound states many of the consequences 
of unitarity still hold true although we do not have any formal proof. This issue should be 
understood more carefully. Infact we find numerically the solution of the SD equation which 
is consistent with the requirements of unitarity is unique and stable. While the other possible 
solutions are numerically unstable. 

Chiral limit of the theory can be consistently defined only if there is no mass UV renor- 
malisation. For this reason alone chiral limit is a critical point since for any other massive 
case the theory is defined only by mass renormalisation. Inspite of this the infrared be- 
haviour of the quark propagator smoothly extrapolates from small mass to vanishing mass. 
Because of this PCAC(Partial Conservation of Axial Current) relations still hold. We find 
chiral symmetric solution of the SD equation and it is also confining due to the absence of 
a pole but it is expected to be an unstable vaccum. 

The techniques of solving SD equations that we have enunciated here is fairly general. The 
massless particle exchange that is considered explicitly has some algebraic simplifications. 
Essentially eq.(J7|) and the trick of eq. ljlUJl can be adopted for any generic relativistic field 
theory to get equations of the type (fTTjl and (JT2J). 

All our considerations above are equally applicable to techincolour models [10J. Indeed 
the existence of many mass scales m, rh, and M is the scenario suggested in the literature, 
here we have a concrete compatible model where it is realised. 
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